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Abstract 

Starting from the Bethe ansatz solution we derive a set of coupled non-linear 
integral equations for the fully packed double loop model (FPL 2 ) on the square 
lattice. As an application we find exact expressions for the central charge and for 
the scaling dimension corresponding to the simplest charge excitation. We study 
numerically the low-lying excitations corresponding to more general perturbations 
of the ground state and discover that the corresponding scaling dimensions are well 
described by the Cartan matrix of s£^. 
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1 Introduction 



The fully packed double loop model (FPL 2 ) is a two dimensional statistical model con- 
structed by decorating the square lattice with two species of loops (black and grey), under 
the constraint that every lattice edge is covered by only one loop and every vertex is visited 
by both types of loops. 

It exhibits a rich phase diagram and provides a representation for previously studied 
models: Ice model [T3], four-colouring model [HJ EEJJ U21 12] , dimer loop model [2H], double 
Hamiltonian walk ^H] and compact polymers [2B]- A generalization of the FPL 2 model, 
which offers a unifying picture of the compact, dense and dilute phase of polymers, was 
proposed in [S] by relaxing the full packing constraint. Another possible generalization 
corresponds to the Flory model of polymer melting [Hl[7j. 

Jacobsen and Kondev mapped the FPL 2 model onto a height model and postulate 
that the long wavelength behaviour is correctly described by a Liouville field theory. 
Making use of the Coulomb gas method they succeed in computing the central charge 
and the critical exponents |U E|. The study of the relation between the Coulomb gas 
representation and the underlying conformal field theory has been initiated in [2*T| . 

Recently we mapped the FPL 2 model onto a 24-vertex model and diagonalized the 
corresponding transfer matrix by means of coordinate Bethe ansatz (BA) PjJ The 
solution holds only when the two loop fugacities are the same. In particular we calculated 
analytically the configurational entropy for the double Hamiltonian walk and the four- 
coloring model. We also confirmed a conjecture on the average length of loops |S]. In 
the particular limit in which the fugacities of the loops are both equal to two, the BA 
equations become rational [2] and look very much like those derived in (22] for mixed 
SU(iV) vertex model for the case N = 4. 

When we made this paper ready for publication, Jacobsen and Zinn- Justin j3j succeed 
in finding a Yang-Baxter structure, for the FPL 2 model, in terms of the affine quantum 
group U q (s£(4)). They reobtain the BA equations following the algebraic approach and 
write exact expressions for the central charge and for the scaling dimensions of the low- 
lying excitations. 

These quantities can be calculated from the leading finite-size corrections to the ground 
[T%1 ITU] and the excited [TBJ [T7j energies. Analytical calculations, for the low-lying ex- 
citations of models solvable by the Bethe ansatz equations, were done in [221 12S]- An 
analytical method, based on the construction of a non-linear integral equation (NLIE) 
for the counting function, which allows the calculation of the central charge, has been 
developed in [33] E2]- This method has been exploited in j33J E01 EU HH3 HUB EH] iu order 
to include excitations. The case of multicomponent BA equations has been investigated 
in 

Here, following [33] , we write a non-linear integral equation for the FPL 2 model in the 
presence of the seam. As an application we compute analytically the value of the central 
charge and of the scaling dimension corresponding to one grey and one black string prop- 
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agating between two points of the lattice. Our results are in agreement with the results 
found in 0] via the field theory. We also studied numerically the low-lying excitations 
generated by varying the number of particles (charge excitations), allowing transi- 
tions from one side to the other of the Fermi seas (umklapp processes) and perturbing 
the surfaces of the Fermi seas (particle-hole excitations) [22] • It turns out that the scaling 
dimensions are well described by a compact formula [23] through the Cartan matrix of 
s£ 4 . Our results coincide with those found in [3]. 

The present work is organized as follows. In section |2] we introduce the FPL 2 model, 
summarize the main results of [I] and define the counting functions. For completeness in 
section El we show in detail how the finite sums entering in the definition of the counting 
functions can be expressed in terms of integrals involving the counting functions. This 
procedure is well known in the literature. See for instance [3*1*1 I34j . In section 0] we derive 
a set of coupled non-linear integral equations for the counting functions, in the presence 
of the seam. An integral expression for the finite-size free energy is derived in section 
The central charge is calculated in section [("] and the scaling dimension associated with 
one grey and one black string derived in section [7| In section |H] we study numerically the 
low-lying excitations. 



The fully packed double loop model (FPL 2 ) is a two dimensional statistical model con- 
structed decorating the square lattice with two species of loops. Each bond of the square 
lattice is covered by only one loop and every vertex is visited by both types of loops. Loops 
of the same type are not allowed to cross. Representing the two species by black and grey 
segments respectively, the fully-packing constraint forces each vertex of the square lattice 
to have one of the six appearances depicted in figure ^ We place the model on a cylinder 
of circumference L and height M with periodic boundary conditions also in the vertical 



packing constraint. Each vertex is visited by both types of loops. The two rightmost 
vertices permit the two species to cross. If we assign an orientation to the loops we obtain 
24 distinct vertices. 



2 The model 




Figure 1: The six vertex configurations of the FPL 2 model that are allowed by the fully- 



3 



direction. The partition function is defined by: 



Z(L,M) = J2< b K 3m t Ib < 9 C 1 ) 



where the sum runs over all the allowed configurations. The exponents Nb and N g are 
the respective numbers of black and grey contractable loops. Since the model is placed 
on a cylinder loops are allowed to wind. We denote by Mb and M g the number of black 
and grey winding loops. We assign the fugacities n& )9 for contractable loops and m bg for 
uncontractable ones. 

In we map the loop model into a 24- vertex model and construct the transfer matrix 
T. The partition function is the sum of the eigenvalues Aj of the transfer matrix, raised 
to the power M: 

Z(L,M) = Tr T M (L) = J^Af (L). (2) 

In the limit in which the height of the cylinder extend to infinite M — > oo the finite-size 
free energy density f(L) can be expressed in terms of the largest eigenvalue Ao(L) of the 
transfer matrix: 

/(L) = ~logAo(L). (3) 

At criticality, the dominant finite-size correction of the free energy density, can be related 
to the central charge c and the bulk free energy /(oo) via the formula [TH| I19j: 

/(L) = /(oo)-^ + o(l/L 2 ). (4) 

The subdominant eigenvalues Aj allow the computation of the scaling dimensions Aj 
thanks to the formula [Tol I17j: 

In P we show that the model is solvable in the particular case in which the fugacities 
are the same (n& = n g and m& = m g ). The transfer matrix is diagonalized using coordinate 
nested Bethe ansatz and an analytical expression for the free energy, in the thermodynamic 
limit, is derived. In particular we find the free energy of the four-coloring model [THIT2"l l2"l 
and the double Hamiltonian walk [TO] and recover the known entropy of the Ice model 
[THj . Let us summarize the main result of pQ: The fugacities entering in the partition sum 
(JTjl are expressed in terms of two phases 9 and a, through the relations: 

n b — n g = 2 cos 7i 9 rrib = m g = 2 cos n a < 9, a < 1 (6) 

and a ^-dependent function is introduced: 

, . sin^(c + xi) 

S c (x,9) = —^ 1 (7) 

sm — (c — xi) 
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With these definitions the Bethe ansatz equations (BAE) assume the form: 



S^uj, 9) L/2 = - Y[Si{wi - uj, 9) Y[ -S 2 {uj - u k , i 

i=l fc=l 

Sl ( Vj , Of' 2 = U S 1 {w i - Vj , 9) J] -Siiyt - v k , t 



i=l 



k=l 

n w 



(9) 
(10) 



exp(-2i7ra) JJSi^ - Uj, 9) Si(wi -v k ,9) -S 2 (w m - w t ,9) = -1. 

j=l k=l m=l 

Every solution of the previous equations, with non coinciding roots, is in correspondence 
with an eigenvalue of the transfer matrix via the expression: 



A 



ex P (i * a) n s l ( Uj , 9f 2 n ^) i/2 n A z - 

m=l I 



3=1 



k=l 



+ ex P (-i ix a) n ft^, #r i/2 n ^ °y i/2 n 



sin^(w m i + 2) 



fc=i 



Sin -7rW m 1 
m=l 2 



In |P the central charge and a set of scaling dimensions were computed analytically 
but not rigorously via the field theory making use of Coulomb gas techniques. In this 
paper we are going to compute the conformal quantities starting from the Bethe ansatz 
solution. These results are included in the independent and recent paper j^j. 

The calculation is based on the formula (£Q) and (0) that relate the central charge and 
the scaling dimensions to the leading finite-size correction of the ground and excited ener- 
gies. An analytical method, based on the perturbation of the Fermi sea, was developed for 
the first time in [221 12H] • Later another method, based on the construction of a non-linear 
integral equation, was developed. Following this method the central charge, of models 
solvable by Bethe ansatz equations, was calculated in [331 IH2] and the scaling dimensions 
in [33 EE EUl EZI- The method has been refined to include string excitations [3T] and 
multicomponents BA equations [231 • Here we are going to follow this last approach. The 
central objects which enter in this machinery are the counting functions Z^ s defined by: 



L 



>l(u,6) + 22 fa ( W ™ ~ u ^)~ /^2fa( u k - U,9) 



m=l 



k=l 



Z LtW (w,9) = } j (j>i(uj-w,6) + } j ^i(vk-w,6) - y^0 2 (w; -w,9) - 2iroi (12) 

3=1 



k=l 



1=1 



Z L>v (v,9) 



+ 2j fa{Wm - v ,Q)-22 fa( Vk ~ V > ^ 
m=l 



k=l 
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where the following function has been defined: 

<j) c {x,6) =i log S c (x,9). (13) 

Notice that the counting functions are constructed by taking the logarithm of the equa- 
tions (JHEII). The key property of Z LjS is that it counts the BA roots, in the sense that 
for each root of the BAE (jSlirUj) we have: 

Z L ,s{^ s ,k) = 27[ h,k k=l,...,n s s = u,w,v (14) 

where I S) & is an integer or half-integer number, depending on the number of roots of each 
type. We will call holes, objects h s ^ that are counted by Z^ s : 

ZL,s{h s , k ) = 2'K\H,s,k k = l,...,n s s = u,w,v (15) 

but are not solution of the system of BA equations (f81ITU|). Where lH,s,k are integers 
or half-integer numbers. The first step towards the exact calculation of the finite-size 
corrections consists in the replacement of the finite sums in the rhs of ()12j) by integrals. 
In the next section we are going to show how this is possible. 



3 From a finite sum to an integral 

We want to find an integral expression for the finite sums entering in the counting functions 
(|12|) for the particular case in which the BA roots A Sj fc are positioned on the real line and 
no real holes are present. These sums have the general form: 

n a 

J^MKk-^d)- (16) 

k=l 

Start by noting that with the help of the counting function it is possible to construct a 
function which is analytic in all the complex plane and has zeros of order one at the BA 
roots. In fact, the following expression satisfies the requirements: 

l + (-l) 5 exp[iZ LiS (^)] (17) 

where we choose 5 = 1 for the case in which the number of roots is such that Z L, s {^s,k) / ^ 
is even and 5 = when it is odd. Applying Cauchy formula we can prove the identity: 

MK k ~ A, 0) = j r d„ M» - A, 9) 1 + M) , exp[iZL>)] (18) 

where T is a close contour which encircles counterclockwise the BA roots X Sy k- A funda- 
mental assumption for the validity of the previous equality is the analyticity of <f) c {z, 0) 
inside and on the contour T. In the following we will use the shorthand notation: 

iZ' Lta (z) (-l) s exp[iZ L , s (z)} 

K S S(Z) = ; -— — . (19) 
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We assumed that the roots of the BA equations were located on the real axis. Therefore 
we choose for T a box centered at the origin, extending from — oo to +00 in the horizontal 
direction and of height 2 rj, and integrate counterclockwise. In the limit 77 — > + only the 
integrals running on the side parallel to the real axis survive and the rhs of (fTHj) becomes: 



du c (/i — A, 9) R s sin) = lim 

77^0+ 



+ 



dx(p c (x - ir/ - X,9) R s j{x - ir?) (20) 



dx 4> c (x + ir/ — A, 9) R s ,s(x + irj) 



The next step consists in the manipulation of the expressions R S; s(x + ir]) and R S) $(x — ir]). 
At a first glance one would say that: 



R St s(x±ir)) = ^log 1 (-1 foxp)Z L J.v± \ti)\ 



[21) 



But this is not always true because the logarithm has a cut on the negative real line. We 
analyze this point in more details focusing on the term R s ,s(x + ir/). Let us study the 
behaviour of the argument of the logarithm in a neighborhood of a solution x of the BA 
equations. To construct the neighborhood take a small real number e and remind that 
77 > 0. The following expansion hold: 



-l^expfi^z + e + i V )} = (77 - ie) Z' L (x) + 0((e + ir]) 2 ). 



(22) 



From this expansion we see that if Z' s L (x) > the argument of the logarithm has a 
positive real part. Thus varying e it will never cross the cut and the equality (|21j) is 
satisfied. On the other hand if y St k is a solution of the BAE for which Z' L s (y s ,k) < (we 
call y S) k a special object) the logarithm will exhibit a discontinuity at e = 0. However the 
expression R s ^(x + ir/) is continuous when x runs from —00 to +00 and so it must be the 
derivative of a continuous function. Thus in the presence of special objects y s ^ relation 
(|2T|) must be replaced with the following equality: 



R s ,s{x + ir/) 



d 

dx 



log 



1 + {-l) s exp[iZ LtS {x + ir])}] + 2vri J] 5(x - y 



(23) 



fc=i 



It can be proved that a similar relation holds for R s> s(x — ir]): 



R 



x-irj) = i Z' LiS (x-ir))+— log [l+(-l)* exp[-iZ L)S {x-irj)]\ -2vri ^ S(x-y Sjk ). (24) 



k=l 



Comparing this equality with (j23j) we see that the special objects contribute with the 
opposite sign and that the extra term Z'{x — ir/) arises. Substituting (J2l?j) and (J2"4~j) into 
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flUJ) yields: 



d/u c (/i — A, 9) R Sj $(jj,) = i / dxcj) c {x — irj — A, 9)Z' L s (x — irj) 



DO 
+ 00 



r+°° d r 

/ dxcj) c (x -\ri - \, 9) — log 1 + (-l) 5 exp[-iZ LiS (x - irj)} 
dx (p c (x + irj-\, °)-^log [l + exp[iZ L , s (x + irj)} 

'•ys n ys 

- 2 tc i ^ ^ciVs^k - irj - K 9) - 2 vr i c (y s ,fc + ir7 - A, 9). 



(25) 



fe=i 



fe=i 



From the definitions (j7ll3H2j) we see that the functions Zl, s and C enjoy the following 
symmetry (since a and 9 are real numbers): 



(j) c (z,9) = <f>c(z,0) = -<f>c(.—z, 
so that we can write the compact formula: 



Z L>s (z) = Z L>s {z) 



dji c (/i — A, 9) R s ^(n) = i / dx <fi c {x — irj — A, 9)Z' L s (x — irj) 



(26) 



(27) 



+00 



d 



— 2ilm / dx<f) c (x + irj — A, 9) — log 1 + (— 1) exp[iZ L s (x + ir?)] 
J— 00 

where we dropped the terms containing the special objects because for the finite-size 
corrections that we will investigate they do not appear. Integrating by parts: 



J dfJL<f> c (fJ, - A, #) R s ,s(fj) = i(pc(x -irj - A, 9)Z L;S (x - irj)\ x x 
dx (p' c (x — irj — A, 9)Zl, s (% — iv) 



x=+oo 
00 



(2f 



— 2 i lm(f) c (x + irj — A, 9) log 1 + (— l) s exp[iZL ;S (x + irj)] 



2ilm / dx(j)' c (x + irj — A, 9) log 1 + (— l) 5 exp[iZ ijS (x + ir/)] 



The line of integration of the first integral can be deformed as long as we stay in the region 
of analyticity of the function C . At the end we get the following integral representation 
for the finite sum (jlfij) : 



"z. r +oc a t 

MKk -\,9) = - — cf>' c (x - A, 8)Z L>s (x) 

k=i J ~°° Zn 



(29) 
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+ lim 2 Im 

+ border terms. 



00 dx 
2^ 



x 



x,e)iog 



■l) s exp[iZ LjS (x + it))] 



In the following, for notational convenience, we will omit the limit symbol and claim 
validity of the equations only in the limit 77 — ► + . 

Notice that in deriving relation (}2*9"|) we made the fundamental assumption that the 
BA roots A s fc were positioned on the real line and that real holes were not present. If real 
holes h S} k are present the equality (J29|) is not valid anymore. Infact the function (|17|) will 
vanish at the real holes h s ^ so that the rhs of (}2*9"j) will contain, beside the contribution 
coming from the BA roots c (A Sj fc — A, 9), also a contribution coming from the BA holes 
Z)<Mft«,k - A, 9). 



4 Non-linear integral equation 

In this section we will derive a set of coupled non-linear integral equations (NLIE) for the 
counting functions. Making use of the key result ()29j) we can reexpress the finite sums 
that enter in the rhs of (|T2*j) in terms of integrals that involve the counting functions: 



Z L , W -2Im$* Q LtV> - 2tt« . (30) 






Ql,u 1 




f 


Q L,w 


l-l 


2 7r a 


Ql,v J 




V 



And we have checked that the border terms vanish. Let us clarify the notation: The 
symbol * means that in performing the matrix product the ordinary product has to be 
replaced with the convolution product. We have defined the derivative of (fT3|) as: 

$ c (*,£)^ J-0 c (x,#) (31) 
and introduced the shorthand notation for the non-linear part: 

Q LiS (x) = log [l + (-if exp[iZ L>s (x + . (32) 
The matrix $ is defined as: 

/ $2 -$1 \ 

$ EE -$j $ 2 -$! . (33) 
V -$! $ 2 / 

In order to solve equation ()3())1 for Zl, s in terms of the non-linear part Ql,s it is convenient 
to work with the Fourier transform: 

= -^(l-l.)- 1 I I -2Im(l-$)- i $ ( Q LtW j (34) 
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_ 15[ _ 

Figure 2: Plot of the terms entering in the second equation of the non-linear system (J38|) . 
Similar plots hold for the first and the third equation. The parameters of the system are 
L = 8, n u = n w = n v = L/2 and 9 = a = 1/4. The curve close to the real axis represents 
the non-linear correction 2 ImG2,s * Ql.s which vanishes in the limit L — > oo. The dashed 
line is the bulk term 2irLz w — 2ira/(l — 9) and the bold line shows the behaviour of the 
counting function Z L w . 





-5{p) (1 - $) _1 (0) [ 2na 



where the Fourier transform and its inverse are defined by: 



F[<j>\{p) = 4>{p) = 
In particular: 



-ipx 



<p(x)dx F- L [(f)](x) = 



1 

2^ 



e ipx 4>(p)dp. 



e- ipx $ c (x } 9)dx 
We will need the following limit: 



sinh(p(c — 1/9)) 
smh(p/9) 



G(p,9) = -(1 - $)^$, lim G(p^) = — — 



1-49 -2 -1 
-2 -AO -2 
-1 -2 1-49 



Working out the terms in ()34|) and computing the inverse Fourier transform 
system of non-linear integral equations (NLIE) for the counting functions: 







( Z u \ 




Zl,w 

\ Zl,v J 


= 2txL 1 




+ lim 2 Im G * 




\z v j 


1 »j-*0+ 





(35) 
(36) 

• (37) 
we get a 

(38) 
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where the bulk quantities are given by: 
z w (x) = — arctan ^tanh(7rx/8) j z u (x) = z v (x) = — arctan (^/2 sinh(7rx/4) j . (39) 

We will need their asymptotic behaviour: 

11 11 

lim z uv (x) ~ -=- exp[— nx/4] lim z w {x) ~ exp[— 7r x/4]. (40) 

x^+oo ' 4 y27T " z^+oo 4 7T 

A plot of the various terms entering in the second equation of ()38|) is shown in figure El 
Similar plots hold for the first and third equation. 

5 Integral expression for the free energy 

In this section we will derive an integral expression for the finite-size free energy. First 
we have to express the free energy as the sum of some function evaluated at the BA roots 
and then applying formula ()29|) convert the finite sum into an integral. For that purpose 
notice that for the distribution of the BA roots corresponding to the largest eigenvalue 
the following equality is satisfied 0: 

flu Tly 

n s m, ^) i/2 n s ^ °y /2 = («) 

3=1 k=l 



which allows to write the finite-size free energy (j3lll|) as: 



HL) = -\ 



log 2 + log cos(2 7r a + s 2 ) + - s\ 



(42) 



where the finite sums s\ and are defined by: 



51 = ^l O g(cos 2 7r^ + sin 2 7r0coth 2 ^|^) (43) 

m=l ^ ' 

52 = arctan (tan7r 9 coth (it 9 w m /2)) . 

m=l 

We will show in the Appendix that only s\ contributes to the value of the central charge. 
In the sequel we will manipulate the sum s\ in order to extract the exact value of the 
central charge. The function entering in the sum S\ (|43|) and its derivative are defined by: 

(p(w) = log (cos 2 n 9 + sin 2 ir 9 coth 2 -^-^-^ ) W(w) = -^-(p(w). (44) 
V 2 / aw 
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Figure 3: Plot of the function H w (x) (dashed line) and of the oscillating non-linear term 
Imlog[l + exp[i Z L)W (x + irj)]]. We refer to a system of size L = 8, phases a = 9 = 1/4 
and parameter rj = 0.01. At the four BA roots Wj the non-linear term exhibits a rapid 
variation. 



At this point remember that in deriving the key result (J29J) which allows the conversion 
of a finite sum into an integral we made the fundamental assumption that the function 
entering in the sum was analytic in the region delimited by the contour integral. And 
moreover that the BA roots where located inside that region. The analyticity requirement 
is not satisfied in this particular case because <p(w) f)44j) is not analytic at the origin. To 
overcome the problem we apply the formalism of sectionEltaking as contour line two boxes, 
one positioned on the right and the other on the left side of the origin, that contain all 
the Wj roots. This is illustrated in figure IU(b). Now, besides the border terms at infinity, 
there will be also border terms at the origin. They are studied in the Appendix. The 
finite sum Si in integral form becomes: 



si 



f dx , , 

= -Vp / — W{x)Z L>w {x) (45) 

J — oo 2lT 

r+°° dx r 
+ 2ImVp / — W{x) log 1 + exp[iZ Lw (x + ir))} 

./-oo 2w I 



+ border terms 

where Vp denotes the principal value. The first term on the rhs of (J45|) can be rewritten 
substituting the expression for the counting function given by the non-linear integral 
equation (|38|). In order to perform the calculations it is convenient to apply Parseval 
formula so that W(x), Z^ w [x) and the non-linear term can be replaced by their respective 
Fourier transform, and then make use of the expression ()34|) for the counting function. 
Performing the calculations in the Fourier space and then applying again Parseval formula 
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via the inverse Fourier transform we get the following integral expression for the sum Si. 



si 



a 



+00 



-LVp / dxW(x) z w (x) + -Vp / dxW(x) 

J -00 1 ~ " J -00 

f +OQ dx r 

+2 Im Vp / —H u (x) log 1 + exp[i Z L , u (x + ir/)] 

Z7T L 



(46) 



00 
+00 



dx 



-2 Im Vp / — H w {x) log 1 + exp[i Z LjW (x + ir])} 

Z7T 



-00 
+00 



dx 



+2 Im Vp / — H v (x) log 1 + exp [i Z L>V (x + ir/)] 
J-00 ^ 

+border terms 



where the functions H s (x) are given by: 

7r sinh(7rx/4) 



H u (x) = H v (x) 



v/2cosh(7rx/2) 



HJx) 



71 



csch(7rx/4). 



(47) 



It is remarkable to notice that the functions H s (x) do not depend on the parameter 9. The 
9 dependence is carried by the logarithmic factors which contain the counting functions. 
Notice that the rhs of ()46j) gives an exact expression for the finite-size sum s\ in terms of 
the finite-size counting functions Zl }S . 

Let us examine the various terms that appear in rhs of ()46|) . The first term is the 
dominant one and is responsible for the bulk free energy. The second one comes from the 
constant part in the NLIE (|38|). It gives a contribution to the finite-size free energy density 
which scales like 1/L. Its contribution exactly cancels with a same contribution coming 
from the border terms. The last three integrals, which exhibit a non-linear dependence on 
the counting functions, will generate the higher order corrections. A plot of the function 
H w (x) and of the non-linear part Im log[l + exp[i Zl, w {x + i rj)]] is shown in figure El The 
non- linear part exhibits rapid jumps at the BA roots Wj. Due to the rapid oscillations, it 
is responsible for the vanishing of the integral in the limit L — > 00. In the next section we 
are going to extract, from these integrals, the corrections of order 1/L 2 which determine 
the value of the central charge. 



6 The central charge 

In order to compute the central charge we have to manipulate the NLIE (J38|) and the 
expressions (I42J) and (j4^|) for the free energy. We make the following change of variable 

x = 4(/i + logL)/7r (48) 
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which depends on the size of the system, and introduce the functions: 

F L>B (x)=exp[iZ Lta (x)] (m 
F +iS (^) = lim L _ +00 F L , s (4(/i + logL)/7r) Q+M = 21m log[l + F +> .(/z + it?)]. 1 ; 

Making use of the asymptotics ()4(JJ) we see that in the limit L — > +00 the NLIE 
reduces to: 











ilog 




,u>(/i) 






^4 








(50) 



+ 21m 



d// G(// — /i) log 



l + F +! „(/i' + iry) 
1 + F +iW (n' + irj) 
l+F +tV {fi' + irf) 



This is a non-linear integral equation for the unknown functions F + S (fi). The dependence 
on the size of the system L has dropped. In particular comparing with (|H8jl we see that 
the terms of order L has been reduced to exponentials. We show now that performing the 
change of variable (|48j) in the expression for the free energy it is possible to extract the 
correction of order 1/L 2 . First notice that the functions H s (x) that enter in the integrals 
have the following asymptotic behaviour: 



7T 

lim H uv (x) ~ = exp[-7rx/4] 

In the new variable \x the leading terms are: 



7T _ 

lim H u>v (A(fj,+ logL)/7r) ~ ^-e" 

so that the expression for the free energy 



lim H w (x) ~ 7T exp[— ir x/4]. 



lim H W (A(^ + log L)/n) 

L— >+oo 



reduces to: 



7T 
I 



(51) 



(52) 



/ (L) ~ / (00) + 



2vrL 2 



+ o(l/^ 2 ) 



(53) 



where the constant C + is defined in terms of integrals of the functions _F +jS (/i): 

/+00 
d/i V2e-» log[l + F +iU (n + i v )} 
-00 
00 

d^2e-" log[l + F +)t0 (/i + i77)] 

/+00 
dfiV2e~^ log[l + F +iV (/x + i77)]. 
-00 



+ 21m 



00 
00 



(54) 
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Notice that the effect of the change of variable (|48|) on the last three integrals of (J46j) is 
twofold: It shifts the oscillating terms so that only one of the two tails of the functions H s 
will contributes to the integral, and in the limit L —>■ +00 it extracts from the function 
H s (A(fi + logL)/7r) the leading term e _At /L. In order to compute the contribution to the 
finite-size corrections coming from the other tail we have to perform the change of variable 
x = 4(fi — log L)/ V. We denote the two contributions by C + and C_. 

Fortunately, it is not necessary to solve the non-linear integral equation (J5U|) in order 
to calculate the constant C + . In fact we can make use of the following formula 



C 4 



-2ReV [ —\og(l + u) 

r,s ~ 

+\^2q+a-° o )q+a-° o ) J 



(55) 



dx G r 



00 
+00 



dx G r 



x 



\x 



where T s is an arbitrary line of integration in the complex plane connecting F +>s (— 00+i + ) 
to -F +s (+oo + i0 + ) and avoiding the logarithmic cut (—00, —1]. This is the multicom- 
ponent generalization of the lemma proved in |3~T] . For completeness we reproduce the 
proof in the Appendix. The limit quantities that enter in the previous formula have the 
following values: 



F +iS (-oo + i 0+) = Q +)S (-oo) =0 s 

F + ' s (+oo + i 0+) = 1 Q+A+°°) = s 

F+ tW (+oo + i + ) = exp[— 2 i n a] Q +jW (+oo) = —2na. 



u, w, v 

u,v (56) 



The integrals involving G. r s can be calculated with the help of the Fourier transform (p?Tj) . 
Summing all the contributions, C + becomes: 



C 4 



r 7T 2 - 


+ 

u 


r 7T 2 


4vr 2 a 2 i 


+ 

w 


r 7T 2 - 


4vr 2 « 2 


. 6 . 


i 6 


h 2 ] 


. 6 . 


, + 2 



e 







71 

~6 



12 a 2 
1-6 



(57) 



The three square brackets refer to the dilogarithm part in formula (|55j) for the three 
families of BA roots, respectively u, w and v. The dilogarithms are worked out using the 
following formula: 

r2 , ,2 



2 Re 



du it" ur 

— log(l + u ) = — - + — 

U D 2 



The same result holds for C_ so that the central charge 

12 a 2 



is given by: 



3- 







(5f 



(59) 
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Figure 4: Arrangement of the roots (black circles) and of the holes (grey circles) in 
the complex plane, a) The u-family. b) The w-family. The parameters are: L — 8, 
n v = n w = 8, n u = 7, a = 9 —>■ 0. The u-family has two real holes. The contour lines 
which encircles the Bethe ansatz roots are shown, a ) Notice that for the w-family the 
contour line encircles also the two real holes, b ) In deriving an integral expression for 
the free energy we avoid the origin splitting the contour line into two boxes. 



which is in agreement with the result found in [U |3| . In particular we see that for the 
four-coloring problem jH E] > which corresponds to the choice 9 = a = the central 
charge is c = 3. For the double Hamiltonian walk {6 = a = 1/2) jHEIl the central charge 
becomes negative c = —3. In general, the shift of the value of the central charge due to 
the presence of the seam has been investigated in [T%1 120j . 



7 The excitations 

The largest eigenvalue of the transfer matrix lies in the sector = L/2 and 

corresponds to roots distributed on the real axis in such a way that the close packing 
constraint is satisfied [Tj: 

^L, s (A Sjfe+ i) - Z LyS (X S)k ) = 2 7r k = 1, . . . ,n s - 1 s = u,w,v (60) 

Subdominant eigenvalues of the transfer matrix can be constructed by perturbing this 
distribution [22] • One way of doing that (particle-hole excitations) is to insert real holes 
h s ^ that will break the close packing of the real roots. Another family of excitations 
(charge excitations) can be constructed by varying the parameters n u , n w and n v that enter 
in the BAE (IHHTUj). A third family of excitations (umklapp processes) can be constructed 
shifting the ground state distribution of the roots. To obtain this last excitation we have 
to replace in (fTlj) the integer numbers l St k, that characterize the distribution of the roots 
corresponding to the largest eigenvalue, by — > l Si k ± 1 for k — 1, . . . , n s . 

In the following we are going to compute analytically, with the help of formula ©, 
the scaling dimension Aj associated with the excitation in which n u is decreased by one. 
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We checked numerically that for the distribution of the BA roots associated with this 
particular excitation relation (|41|) is satisfied so that the finite-size energy associated with 
the excited eigenvalue Aj(L) is given by: 



log Aj (L) = log 2 + log cos(2 tc a + s 2 ) 



1 

- Si 
2 



(61) 



where the sums Si and s 2 are defined by ()43|) . We show in the Appendix that also in 
this case the contribution to the value of the scaling dimension comes exclusively from 
the sum s\. Thus, the machinery developed in section El for the calculation of the central 
charge can be applied, with small changes, also in this case. Solving numerically the BA 
equations for the excited state we see that two real holes appear for the w-family. This 
is shown in figure 0] (a). There are no real holes for the other two families of roots. To 
denote the position of the real w-holes we introduce the positive real numbers h% for the 
hole located on the positive real axis and h 2 for the hole on the negative axis. Also in 
this case, following the steps of section 0J we can write a non-linear integral equation for 
the counting function. For that notice that the sums, entering in the counting functions 
(112)1 . run over the BA roots. On the other hand the contour integral is a box centered 
around the origin and extending from — oo to +oo, so that it will encircle also the two 
holes located at the border of the w-roots (see figure E] (a)). At the end their contribution 
must be subtracted. The analog of Eq. (jSOjl becomes: 

Z L ,u( X ) 

Z L)W {x) 

Zl,v( X ) 



+ 




(62) 



/+oo 
dx' G(x' — x) log 
-oo 



1 - exp[iZ L>u (x' + ir))} 
1 + exp[iZ L , w (x' + irj)} 
1 + exp[iZ L , v (x' 



Notice that in the non- linear term of (|52j) there is a minus sign in front of exp[i Zl u ]. This 
is because for the excited state we are investigating the counting function Zl, u evaluated 
at the BA roots and holes is an even multiple of ir. With this change of sign it turns 
out that the function (1 — exp[i Zl, u ]) has a zero of order one at the BA roots and holes, 
a necessary condition to apply the machinery developed in section El Repeating the 
calculations of section |U we can write the analogue of (|5Sj) for the particular excited state 
we are investigating: 



Z L ,u{x) 
Z l ,w{x) 

Z L ,v( X ) 



2ttL 





S u (x + h 2 ) + S u (x - hi) 
S w (x + h 2 ) + S w (x - hi) 
S v (x + h 2 ) + S v (x - hi) 
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+2Im / dx' G(x' - x)\og 



1 - exp[iZ L)U (x' + ir))} 
1 + exp[iZ L)W (x' + ir})} 
1 + exp[iZ L>v (x' + ir])] 



(63) 



where the vector S can be expressed in terms of the inverse Fourier transform of the 
matrix G(p, 0) (J37|): 




+oo 



dp 



smxp 
p 




(64) 



In particular it turns out that the derivative of S is related to G: 

SJx) 



-f- I S w {x) 
ur 1 SJx) 



2tt 



G n (x) 
G 2 i(x) 
G 31 (x) 



(65) 



From expression and (|H7|) we read the limit values of S: 



(l-40)/(4-40) 
±7T I l/(-2 + 2 0) 
l/(-4 + 40) 

and taking the limit L 



7r a 




Performing the change of variable 
reduces to: 



(66) 

+oo, equation (J53)l 





" F + 


,M " 




ilog 


F + 




-ay 






M _ 





(1-0) 



(67) 



(1 -40)/(4-40) 
- 77 | l/(-2 + 2 0) 
l/(_4 + 40) 

/+oo 
dfi' G(// - n) log 
-oo 



S v (fi - fi hl ) 

1 - F +iM (/i' + ir]) 
1 + F+^lfjf + iri) 
1 + F +tV (fj,' + i V ) 



which is the analogue of (JHOJ) for the excited state. We will need the following limit values: 

(68) 



F +iS (-oo + iO) =0 
F + , u (+oo) =exp[i7r(l-2 0)] 
F +)UI (+oo) = exp[i7r( — 1 — 2a 
F+J+oa) = 1 



Q +)S (-oo) 
Q+,J+oo) 
Q +!W (+oo) 
Q +t J+oo) 



s 
-2tt0 

= 7r(— 1 - 2a + 
0. 



U, W, V 
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The excitation we are studying modifies also the integral expression for the finite-size free 
energy. In particular in the integral expression for the sum s x (|46|) the following extra 
term will appear: 

y °° dx 

—W(x) (S w (x + h 2 ) + S w (x - h)). (69) 

In the following we will compute the contribution to the finite-size correction coming from 
the hole located at position hi. Similar calculations hold for the hole positioned at — h 2 . 
Working out the integral we obtain: 

+ -^ w(x) sj x - hl) ^J^^l±tl!4) (70 ) 

oo 2tt kj wK 11 6 ^cosh(ft! ir/4)-l/V2' 
which has the following asymptotic behaviour: 



Um ^cosh( ftl ,/4) + l/^\ 
6 \ v cosh(/i 1 vr/4) - l/y/2 J 



xvr/4]. (71) 

Following the steps of section El we see that, in the limit L — > +oo, the finite-size energy 
of the excited state reduces to: 

log Ai(L) = bulk + C \ JrC ~ + o(l/ L) (72) 

Z 71 L 

where the constant C + is now given by: 

C + = -27rV^exp[-^ Al ] (73) 

/+oo 
dnV2e~^ log[l - F +iU (/i + i 77)] 
-00 

/+00 
dn2e~ IM log[l + F +tW (/i + iT])} 
-00 

/+00 _ 
dfi V2e-» log[l + F +iV} (ji + ir})]. 
-00 

The bulk term in (|T2*|) coincides with the bulk term of the ground state eigenvalue (p)3*|) . 
The dominant finite-size correction has a different structure. In fact comparing (|T3*|) with 
(p>4"j) we see that there is an extra contribution — 2n \/2exp[— /i/jj due to the presence of 
the hole positioned at and there is a minus sign in front of F +tU . A similar expression 
can be derived for C_ performing the change of variable x = 4(/z — logL)/7r which selects 
the hole located at — fih 2 - Now we are going to rewrite (|T3*|) in a fashion that allows the 
application of formula (jSHJ)- For that notice that the hole located at solves the first 
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equation of the system ()67|) for some integer number l hl : 

7r a ir(l — 4 1 



2vrl 



hi 



■V2exp[-fi hl ] 



S u (0) 



(1-6) 4(1-6) 
— ^-Sui/j, - fi hl ) log[l - F +jU (ii + i 77)] 

2lm [ 7r--rS w (n - fi hl ) log[l + F +:W (/j, + i 77)] 

21m / — ^(M-MfcJlogtl + ^^ + i^)] 

J -00 2 7rd/z 



(74) 



where according to (JMj) S^O) = 0. Now, substituting the expression for the exponential 
given by relation (j74j) into the rhs of (J73J) we get the following expression for C + : 

2n 2 a 7i 2 (l-46) 



C 4 



4vr 2 I, 



- 21m 

- 21m 

- 21m 



-00 
+00 



■00 
+00 



1-6 2(1-6) 

r d 
d\x \f2e~ 11 + — S u (n 
L d/i 



(75) 



log[l-F +ltt (/i + i77)] 



dfj, 



d 



2e~" + -^S w (fi - /^JJ log[l + F+^ifi + in)] 
V2e~^ + -£-S v (ji - n hl )] log[l + F +iV (/i + I77)]. 



Notice that the expression in square brackets is the derivative of the first four terms on 
the rhs of (joTj) so that we are in the position to apply formula (|55p. Showing all the terms: 

C+ = 4irl hl + - — — — (76) 



7T~ 



1 

L3 



1-6 2(1 
'1-2 0) 2 



1 

6 



■l-2a + 





- r 




+ 

- w 


-~ 6. 


+ q) 

V / 



For C_ we obtain the same expression with 1^ replaced by I^ 2 . The three square brackets 
refer to the dilogarithm part in formula ()55|) for the three family of BA roots, respectively 
u, w and v. The dilogarithms are worked out using formula (J55j) and the limit values 
f!68|) . In order to compute the contribution coming from the w-term we make use of the 
formula: 

du , . 7T 2 

-2Re/ —log(l-u) = _ + _-7ra;. 
Jo u 3 2 

The term g in the expression (J75J) gives the contribution coming from the non-dilogarithm 
part of formula (|55|l . Working out we find: 

0(3 + 4 a 2 + 8 a - 4 a - 5 6 + 5 fl 2 ) 



2(1 
20 



(77) 



We have all the ingredients to apply formula (JHJ). Substituting the values of C_ and C + 
in (|72j) and reminding that the central charge is given by (J59j) after some algebra we get 
simple expressions for the scaling dimensions: 

A + = -l hl + ^- A_ = -I fe + i^ A = A+ + A_ = ^ (78) 

where the a dependence has dropped. From numerical calculation we see that the exci- 
tation under investigation corresponds to the choice = lh 2 =0. 

From this result we can compute the scaling dimension corresponding to one black and 
one grey string propagating between two points on the lattice jl]. Call Gi(r) the two string 
correlation function associated with defects configurations where a single oriented black 
loop segment and a single oriented grey loop segment propagate from a vertex located at 
the origin to a vertex positioned r rows above it. If we denote by <j) the operator which 
creates the defect, the correlation function, in the transfer matrix formalism is given by: 

d(r) = (0(r)0(O))= lim llVT^0(r)T!l 1AO 0(O) (79) 

N^OO Zj 

= ^e- ( ^ (L) - Eo(L))r (O|0(r)|n)(n|0(O)|O) 

n 

where T 0> o,o is the block of the transfer matrix corresponding to the ground state sector 
and T_ 10i o is the block corresponding to the excited sector in which the number of u 
particles is decreased by one. Notice that the sum runs over the eigenvalues e~ E ^ L ) of 
T-1,0,0 an d e~ E °( L ^ is the largest eigenvalue of the transfer matrix. In the limit r ^> 1 the 
correlation exhibits, in the cylindric geometry, an exponential decay: 

Gi(r) ~ e - {El - Eo)r (80) 

where e~ El ^ L > is the largest eigenvalue of the matrix T_ lj0)0 - In the bulk (L — > oo) the 
correlator exhibits a power law behaviour: 

Gi(r) 



If the model is conformaly invariant the energy gap and the critical exponent can be 
related via the formula (see also formula (JHJ)): 

Xl= Km (E\{L) — E (L)). (82) 

L— >oo Z 71 

It is important to notice that since generates the configurations in which a black 
and a gray string run from the first to the rth row of the cylinder, black and grey winding 
loops are forbidden inside this region so that there is no need to introduce the seam in order 
to assign them the right fugacity |2Z] • Moreover an eventual seam would associate 
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a) 



b) 
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1 — rtr 
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Figure 5: Arrangement of the roots (black circles) and the holes (grey circle) in the 
complex plane for the w-family. a) Arrangement corresponding to the ground state. The 
parameters are: L = 8, n u = n v = n w = 4, a = 9 = 0. Notice that the ground 
state possesses complex holes, b) Arrangement corresponding to the excitation labelled 
by m = (—1,0,0) and e = (1,0,0). Notice that the charge excitation m = (—1,0,0) 
generates two real holes at the border of the Fermi sea (see figure 0] (a)). The umklapp 
process e = (1, 0, 0) shifts all the root distribution to the right. 



a phase, to loops spiralling around the cylinder, proportional to the winding number. 
Performing the calculation that lead to (|7Hj) with a = 9 for the largest eigenvalue and 
a = for the subdominant one we get for the exponent governing the two string correlation 
function: 

1-9 9 2 

x ^ — -—e (83) 

in agreement with the result found in [3] via the field theory. Notice that the second term 
in the rhs of (|83j) is 1/12 of the shift of the central charge (see the formula (|59|) for the 
central charge). 



8 Numerics 

In the previous section we have shown in detail how to compute analytically the scaling 
dimension associated with the excitation in which the parameter n u is decreased by one. In 
general, an excitation will be characterized by a particular arrangement of holes and roots 
in the complex plane. If the arrangement is known the corresponding scaling dimension 
can be calculated going through the machinery that we have described. In this section 
we are going to study numerically a larger set of excitations. 

To classify the excitations in a compact manner we introduce three vectors e = 
(e u ,e w ,e v ), m = (m u ,m w ,m v ), I ± = (l u ,l w ,l v ). The components of m are integer 
numbers which count the change of the "number of particles" in the excited state with 
respect to the ground state. The components of e count the total number of transitions 
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of the "particles" from one side to the other of the three Fermi seas. The non-negative 
integers If describe excitations in the vicinity of the Fermi surfaces. They count the num- 
ber of holes that are generated by moving the roots, positioned at the Fermi surfaces, to 
the next vacant positions. As example, figure |3] shows the arrangement of holes and roots 
corresponding to the excitation labeled by m = (—1, 0, 0), e = (1, 0, 0) and P* 1 = (0, 0, 0) 
for a system of size L = 8. 

Table Q shows a set of charge excitations labelled by the vector m. The scaling 
dimensions have been estimated by solving numerically the BA equations for systems of 
increasing size and making use of formula (jSJ). Notice that in applying formula (0) the 
largest and the subdominant eigenvalues were computed keeping the value of a constant 
and not varying it as in the derivation of (|83|) . From the data we see that the scaling 
dimensions do not exhibit a dependence on a. The data for m = (—1,0,-1) and m = 
(— 1, 0, +1) merit a special comment. Solving numerically the BA equations for a generic 
value of 6 and setting a = we see that two particles at the border of the w-family go 
to infinity. Specifying the value wi = ±oo in (fTUj) and observing that for the function (|7J) 
lima.^-1-oo S c (x,8) = — exp [=pi tt 9 c] it follows that (jlOj) is satisfied only for a = 0. Thus 
the excitations m = (—1, 0, —1) and m = (—1, 0, +1) we have founded only for a = 0. 

Table |3] shows the scaling dimensions for a set of excitations constructed perturbing 
the ground state by a combination of charge excitations labelled by m with umklapp 
processes labelled by e. We have checked that for this set of excitations in which the 
umklapp is involved the scaling dimensions exhibit a dependence on a. 

Table |2] shows the scaling dimensions for a set of excitations constructed combining 
charge excitations m with particle- hole excitations labelled by I 1 * 1 . In order to construct 
the particle-hole excitations we insert a hole at the very end of the Fermi sea. Notice 
that there are three Fermi seas, one for each type of particles, and each Fermi sea has two 
surfaces that can be perturbed independently. We distinguish among the two surfaces 
using an upper index for the vector I ± . From the data we see that the scaling dimensions 
for this class of excitations do not depend on the phase a. Notice that the scaling di- 
mension for the excitation labelled by m = (—1, 0, 0) and I + = (1, 0, 0) can be calculated 
analytically setting in formula (J7HJ) I/ ll = — 1 and lh 2 = 

The set of excitations shown in table ^ and H3 does not cover all the possibilities. In 
fact we did not succeed in finding the excitation labelled by e = (+1, +1, 0) and the one 
labelled by m = (+1, —1,0). The numerical data, of the low-lying excitations we could 
construct, suggest that the scaling dimensions A, are described by the formula: 

A(m, e,I ± ,e)~m i Cm + r ^ e* C" 1 (e - e ) + £ 1+ + £ lj (84) 

s s 

where the background charge is given by eo = (0, 2 a, 0), and C is the Cartan matrix of 
sU / 2 -1 

-1 2 -1 | (85) 
0-12 
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1X1 







L = 100 


L = 200 


L = 300 


L — oo 


(-1,0,0) 
(+1, U, UJ 


0.25 





0.375018 


0.375004 


0.375002 


0.375 


[—1, u, UJ 
(+1,0,0) 


0.25 


0.25 


0.374946 


0.374986 


0.374994 


0.375 


(0,-1,0) 

(-1,-1,0) 


0.25 


0.25 


0.375018 


0.375005 


0.375002 


0.375 


(-1,0,-1) 

(-1,0, +1) 


0.25 





0.750368 


0.750092 


0.750041 


0.75 


(-2,0,0) 
(0, -2,0) 


0.25 





1.499419 


1.499854 


1.499935 


1.5 


(0,-2,0) 


0.25 


0.25 


1.499421 


1.499854 


1.499935 


1.5 



Table 1: A list of charge excitations labelled by the vector m. We group in the same 
box excitations that have exactly the same excited eigenvalue. The scaling dimensions 
are estimated by solving numerically the BA equations for systems of increasing size L 
and making use of formula Notice that the value of the scaling dimension does not 
depend on a. The last column shows the dimension predicted by formula ()84j) . 



m, I + 


9 


a 


L = 100 


L = 200 


L = 300 


L = oo 


m= (-1,0,0) 
I+ = (1,0,0) 


0.25 





1.373801 


1.374700 


1.374867 


1.375 


m=(-2,-2,-2) 
I+ = (1,0,0) 


0.25 





2.495597 


2.498891 


2.499506 


2.5 


m=(-2,-2,-2) 
I+ = (1,0,0) 


0.25 


0.25 


2.4966007 


2.499141 


2.499617 


2.5 


m= (-2,-3,-2) 
I+ = (0,1,0) 


0.25 





2.871114 


2.874015 


2.874560 


2.875 


m=(-2,-3,-2) 
I+ = (0,1,0) 


0.25 


0.25 


2.873027 


2.874490 


2.874771 


2.875 



Table 2: A list of excitations constructed combining the charged excitation labelled by 
the vector m with the particle-hole excitation labeled by I + . The scaling dimensions are 
estimated by solving numerically the BA equations for systems of increasing size L and 
making use of formula ©• In order to create the particle-hole excitation we insert a hole 
between the two particles that are located at the very end of the Fermi sea. Notice that 
for this class of excitations the scaling dimensions do not depend on the phase a. The 
last column shows the dimensions predicted by formula ()84jl . 
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(X 


L=100 


L=200 


L=300 


L — OQ 


e= (-1,0,0) 
m — (— 1, U, UJ 


0.25 





1.373677 


1.374668 


1.374852 


1.375 


e=(-2,0,0) 
m — (—4, -4, — 4J 


0.25 





9.922879 


9.980360 


9.991231 


10 


e=(-l,0,+l) 
m - {— 2, -2, -2) 


0.25 





2.829248 


2.832285 


2.832863 


2.83 


e=(+l,0,+l) 

m — [—2, —2, —2) 


0.25 





4.154368 


4.163582 


4.165295 


4.16 


e=(0,+l,0) 
m= (-1,-2,-1) 


0.25 





2.087041 


2.084259 


2.083745 


2.083 


e=(-l,+l,0) 
m= (-2,-2,-2) 


0.25 





2.875721 


2.875168 


2.875073 


2.875 


e= (+1,0,0) 
m = (-2, -2, -2) 


0.25 


0.3 


2.097194 


2.099292 


2.099684 


2.1 


e= (-1,0,0) 
m = (-2, -2, -2) 


0.25 


0.3 


2.893517 


2.898372 


2.899276 


2.9 



Table 3: A list of excitations constructed combining a variation in the number of "par- 
ticles" labelled by m, with transitions, from one side to other of the Fermi seas, labelled 
by e. The scaling dimensions are estimated by solving numerically the BA equations for 
systems of increasing size L and making use of formula (JSJ). Notice that for this set of 
excitations in which e ^ (0, 0, 0) the scaling dimensions exhibit a dependence on a. The 
last column shows the dimensions predicted by formula (j84j) . 



which confirms the results proposed independently in A formula of this type was 
derived for the first time in [23] for the conformal dimensions of models solvable by 
multicomponent Bethe ansatz. A rigorous proof of the fact that the scaling dimensions 
associated with the charge excitation and the umklapp process are essentially inverse of 
each other has been given in [22] • A similar expression for the conformal dimensions, in 
which C is replaced by the Cartan matrix of si^, has been found in [21] for a loop model 
on the hexagonal lattice. 

9 Conclusions 

We have derived a non-linear integral equation for the FPL 2 model on the critical part of 
the solvable line. As an application we have extracted exact expressions for the central 
charge and for the critical exponent associated with one black and one grey string propa- 
gating between two points on the lattice. We have also studied numerically the low-lying 
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excitations that can be constructed perturbing the Fermi seas in a more general manner 
(charge, umklapp and particle-hole excitations) and found that the scaling dimensions can 
be classified in a compact way through the Cartan matrix of sl^. The calculations have 
been done without having the knowledge of the underlying group structure discovered in 
[3] and constitute a confirmation of their results. It would be interesting to continue our 
study by considering strings excitations that can be generated selecting solutions of the 
BA equations in which some of the roots have an imaginary part different from zero. 

A Appendix 

A.l The sums s\ and 

Here we are going to show that the value of the central charge, which gives the dominant 
finite-size correction to the free energy density (J42j) . is determined exclusively by the sum 
s\. The sum s\ and S2 has been defined as 



E. . / 2/1 2 ft 1 2 ts CC \ / \ 

(PliWm) ifiyx) = log ^ cos 7r0 + sin vrfc'coth — — J (86) 



771=1 



S2 = /J ¥>2(w m ) V ? 2(^) = arctan(tan(7r 9) coth(7r 9 x/2)). (87) 

777=1 

The function <fi(x) has a singularity at the origin which, in deriving the integral expression 
for the free energy, will generate border terms. In particular it can be shown that: 

1 r 1 °" 

lim -Im^i(x + ir/)log 1 + exp[iZ L , w (x + i 77) = 2(log2 + logcos(Z L ,«(0)/2)). (88) 

77->0+ 7T L J 0+ 

We checked numerically with fifteen digits of accuracy, that, for the ground state and the 
excited state under investigation, the term depending on S2, in the expression for the free 
energy ()42j) . exactly cancels with the previous border term: 

log cos(2 7T a + s 2 ) - log cos(Z Li „,(0) /2) = (89) 

so that the finite-size corrections to the free energy density are determined exclusively by 
the sum si. 

A. 2 Lemma 

Here, for completeness, we reproduce, following j^I], the prove of formula (|55j). 
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Hypothesis: 

Assume that F(x) satisfies the non- linear integral equation (NLIE): 



— i log F(x) 



(x) + 



dyG{x-y)Q{y) 



where the non-linear term Q(x) is defined by: 

Q(x) = 2 Im log[l + F(x + i rj)] 



(90) 



(91) 



assume also that the functions <p(x) and G(x) are real and that G(x) is symmetric, inte- 
grable and peaked around the origin. 
Thesis: 

Then the following identity holds: 



dx 4>'(x) Q(x) 



-2 Re / — log[l + d 
Jr u 



(92) 



g 2 (+oo) -g 2 (-oo) 



+00 



dxG(x) 



where T is a contour in the complex u-plane from F(— oo + i0 + ) to F(+oo + i0 + ) avoiding 
the logarithmic cut (—00, —1]. 
Proof: 

Deriving (}9T)|) and substituting the expression for <p'(x) into the lhs of (|9*2*|) we get: 



dx <f)'(x) Q(x) 



f+°° F'(x) 
-21m / dxi— -ry log[l + F(x + ir])] 
J -00 r \ x ) 

dyG\x-y)Q{y) Q(x) 



(93) 



dx 



The first term on the rhs of (}9~3*j) originates the dilogarithmic term in the rhs of ([92)1 . Let 
us examine in detail the double integral l(a, b) where < a < b: 



l(a,b) 



dx 



dx 



dyG'{x-y)Q(y) Q(x) 



(94) 



6 

—a 



+b 



dyG'(x-y)Q(y)+ / dy G'(x - y) Q(y) Q(x) 



where we have used the fact that G'(x) is an odd function. Integrating by parts: 



I{a,b) 



dx 



+ 



+a 



dx 



G(x + b) Q{-b) - G{x + a) Q{-a) + / dy G{x - y) Q'{y) Q(x) 

J-b J 

r+b 

G{x-a)Q{a)-G(x-b)Q(b)+ / dy G(x - y) Q\y) Q(x) (95) 
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In the limit lim a ^ +QO [lim fo _^ +00 I(a, b)} the integrals and the terms depending on b, inside 
the square brackets of ([95]). vanishes because lim^-i-oo G(x) = 0. For the terms depending 
on a we can write: 

/+a i r+oo 

dx[TG(x±a)Q(Ta)Q(x)} =T-Q 2 (Too) dxG(x) (96) 
-a 2 

where we have used the fact that G(x) is symmetric and peaked around the origin together 
with the assumption that Q(x) tends to some asymptotic value in the limit x — > ±00. 
The generalization to the multicomponent case jSl] is straightforward. 
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